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ABSTRACT 

A finite element formulation is derived for the scattering and radiation of 
acoustic waves from submerged elastic structures. The formulation uses as fundamen- 
tal unknowns the displacement in the structure and a velocity potential in the fluid 
Symmetric coefficient matrices result. The outer boundary of the fluid region is 
terminated with an approximate local wave-absorbing boundary condition which assumes 
that outgoing waves are locally planar. The finite element model is capable of pre- 
dicting only the near-field acoustic pressures. Far-field sound pressure levels may 
be determined by integrating the surface pressures and velocities over the wet bound 
ary of the structure using the Helmholtz integral. Comparison of finite element 
results with analytic results show excellent agreement. The coupled fluid-structure 
problem may be solved with general purpose finite element codes by using an analogy 
between the equations of elasticity and the wave equation of linear acoustics. 

INTRODUCTION 

There is a variety of practical engineering problems which cannot be addressed 
using only the separate disciplines of structural analysis and acoustics; they must 
instead be treated by formulating a coupled structural-acoustic problem. For exam- 
ple, for aircraft (Ref. 1), automobiles (Ref. 2), and railroad cars (Ref. 3), cou- 
pled analyses are being investigated to understand better the interior noise problem 
so that such noise can be reduced. In aerospace vehicles, the vibrations of fluid- 
filled tanks are of interest (Ref. 4). In aircraft hydraulic systems (Ref. 5) and 
shipboard piping systems (Ref. 6), the dynamic behavior of (and the transmission of 
sound in) fluid-filled piping systems has been analyzed. Other important naval prob 
lems involve the vibration of underwater structures such as rudders and propellers 
(Ref. 7-8), the shock response of submerged structures (Ref. 9-13), and the scatter- 
ring of sound waves from underwater elastic structures (Ref. 14-17). 

The commonality among all these problems is the mathematical model. The struc- 
ture, if it can be assumed to remain elastic, behaves according to the classical 
theory of elasticity (Ref. 18) and the various approximate engineering theories for 
beams, plates, and shells. The fluid is generally treated as an acoustic medium 
(Ref. 19-21), a fluid whose pressure p satisfies the scalar wave equation 


y2p = p/c^ 


( 1 ) 


where c is the speed of sound in the fluid. The boundary condition at a fluid- 
structure interface can be obtained from momentum and continuity considerations: 


3p/9n = -pu 


( 2 ) 
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where n is the normal at the interface, p is the mass density of the fluid, and 
u n is the normal component of fluid particle acceleration. 

Dynamics problems involving the interaction between an elastic structure and an 
acoustic fluid have been formulated for finite element solution (Ref. 22) by using 
either pressure (Ref. 19,23,24) or fluid particle displacement (Ref. 14,25-27) as 
the fundamental unknown in the fluid region. In three dimensions, the pressure and 
displacement formulations result in, respectively, one and three degrees of free- 
dom per finite element mesh point. Thus the pressure approach has the advantage of 
fewer unknowns and a smaller overall matrix profile or bandwidth. On the other hand, 
the displacement approach results in symmetric coefficient matrices (in contrast to 
the pressure formulation, for which the matrices are nonsymmetric) and a fluid-struc- 
ture interface condition which is easier to implement with general purpose finite 
element computer programs. However, the displacement approach also suffers from the 
presence of spurious resonances (Ref. 27), a situation which can be bothersome in 
time-harmonic problems, either forced or unforced. Recently it was shown (Ref. 28) 
that the principal disadvantage of the pressure formulation, nonsymmetric coefficient 
matrices, can be removed merely by reformulating the pressure solution approach so 
that a velocity potential rather than pressure is used as the fundamental unknown 
in the fluid region. For some situations, particularly steady-state problems involv- 
ing damped systems and time-dependent problems, significant computational advantages 
result . 

The principal goal of this paper is to present in detail the symmetric velo- 
city potential formulation for application to the specific problem of acoustic 
scattering from submerged elastic structures. Previously (Ref. 28), the symmetric 
potential formulation was described only in general terms for a wider class of 
fluid-structure interaction problems with no details concerning specific types of 
applications (such as vibrations, shock response, or acoustic scattering). 

The scattering approach described here has advantages over the displacement 
formulations for the reasons already given. Finite element modeling of exterior 
fluid regions also has advantages over the use of approximate theories such as the 
doubly asymptotic approximation (DAA) (Ref. 15). In this case, the primary trade- 
off is between the large, banded matrices which finite element models generate and 
the smaller, densely-populated matrices which the DAA generates. This trade-off 
often favors the finite element approach for long structures like ships which are 
"naturally banded." Similar trade-offs arise if the finite element approaches are 
compared with T-matrix methods (Ref. 16). T-matrix approaches, however, are not 
yet available in general purpose codes capable of handling arbitrary three-dimen- 
sional geometries. 

From an engineering point of view, it is convenient to be able to make use of 
existing general purpose finite element codes (such as NASTRAN, among others), be- 
cause of their wide availability, versatility, reliability, consultative support, 
and abundance of pre- and postprocessors. Thus the next section summarizes an ana- 
logy between the equations of elasticity and the common field equations of classi- 
cal mathematical physics (including the wave equation). This analogy allows the 
coupled structural-acoustic problem to be solved with standard finite element codes. 

Subsequent sections of the paper will develop the formulation of the scatter- 
ing problem for elastic obstacles with and without fluid inside. Examples will be 
shown for both cases. 
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STRUCTURAL ANALOGIES FOR SCALAR FIELD PROBLEMS 


Since we wish to solve the coupled structural-acoustic problem using standard 
finite element codes (which were developed principally for structural analysis), we 
summarize here the application of such codes to various nonstructural field problems 
(Ref. 24,29). 

Many linear problems in mathematical physics involve the solution of an equa- 
tion obtained by specializing the general form 

V 2 <|)+g=a<£+b<|> (3) 


where V 2 is the Laplacian operator; dots denote partial time differentiation; the 
functions g, a, and b are, in general, position-dependent; and the unknown scalar 
function <j) depends on both position and time. 

Special cases of Eq. (3) arise in such diverse applications as heat conduction, 
acoustics, electrical and magnetic potential problems, torsion of prismatic bars, 
potential fluid flow, and seepage through porous media. Several common special 
cases are listed here: 


Laplace’s equation: 

V 2 <|> = 0 

(4) 

Poisson f s equation: 

V 2 + g = 0 

(5) 

wave equation: 

V 2 <t> = '<ji/c 2 

( 6 ) 

heat equation: 

• 

kV^<j) 4- q - pc<J> 

(7) 

telegraph equation: 

• 

3 2 <j>/3x 2 = LC<j> + RC(|> 

( 8 ) 

Helmholtz equation : 

V 2 <|> + k 2 4> = 0 

(9) 


Most boundary conditions likely to be encountered in connection with Eq. (3) 
will probably be special cases of the general form 

* *• 

aj 3<|>/8n + <t> + £3 <|> + (j> + a^ = 0 ( 10 ) 


where n is the outward normal at the boundary. For example, in heat conduction 
problems, a boundary with a prescribed temperature function satisfies the Dirichlet 
condition 


<t> “ ^o 


(ID 
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and a perfectly insulated boundary has the Neumann condition 


3<}>/9n = 0 


(12) 


In free surface flow problems, the linearized free surface condition on the velocity 
potential is (Ref, 30) 


♦ + §o *,z = 0 (13) 

where g Q is the acceleration due to gravity, the free surface is the plane z = cons- 
tant, and commas denote partial differentiation. In one-dimensional radiation prob- 
lems, the plane wave radiation condition that the velocity potential must satisfy 
at a non-reflecting boundary is 


»,n + 4>/ c = 0 


(14) 


where c is the wave speed. 

An example of a boundary condition not of the general form of Eq. (10) is the 
condition (2) which must be satisfied at an accelerating boundary of a fluid. 

According to the analogy (Ref. 29) between Eq. (3) and the Navier equations of 
classical elasticity, Eqs. (3) and (10) can be solved with elastic finite elements 
using the following procedure: 

1. Select one of the three Cartesian components of displacement (or the 
z-component in cylindrical coordinates) to represent the scalar field variable <j>. 
Constrain all other displacement components everywhere in the field. 

2. Model the domain of interest (either 2-D or 3-D) with finite elements 
having material constants satisfying 


E e - a G e , p e - a G e 


(15) 


where a is the variable appearing in Eq. (3), and E e , G e , and p e denote the Young’s 
modulus, shear modulus, and mass density assigned on the material card to the finite 
elements. The subscript "e" has been added to emphasize that these constants are 
merely numbers assigned to the elements and may bear no resemblance to any actual 
material properties associated with a particular application. The dimensionless 
constant a in Eq. (15) should, for 3-D problems, be chosen large enough to make otfl 
numerically indistinguishable from ot. For 2— D problems, a should be small, but not 
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so small that 1+a is numerically indistinguishable from unity. Thus, on most compu- 
ters , 


10-5 

(2-D) 

1020 

(3-D) 


(16) 


The shear modulus G e can be selected arbitarily. The finite elements eligible for 
use in the model are those derived from classical elasticity theory rather than from 
the engineering theories involving beams, plates, or shells. Thus, for 2-D problems, 
the plane stress membrane elements are appropriate. For 3-D problems, the solid 
elements (e.g., the isoparametric or the axisymmetric solids) should be used. 

3. Apply to the unconstrained degree of freedom (DOF) at each grid point in 
the region a "force” given by 


F = G e g V 


(17) 


where V is the volume assigned to the point and g is the function appearing in Eq. 
(3). For problems for which the function g in Eq. (3) is independent of position 
(as, for example, in the classical St. Venant torsion problem), this load may be 
specified conveniently by applying to the "structure" a gravitational field for which 
the acceleration due to gravity g Q satisfies 


Pe So ^e & 


(18) 


4. Connect between ground and the unconstrained DOF at each grid point in 
the region a scalar dashpot whose damping constant (the ratio of damping force to 
velocity) is G e bV, where b is the function appearing in Eq. (3) and V is the volume 
assigned to the point. 

5. Enforce the boundary condition (10) by applying to the unconstrained DOF 
at each grid point on the boundary of the region a "force" given by 


F — — GgA (a£ 4^ F U 3 4* F* a 4 4* a^)/a^, 


a^ * 0 


(19) 


where A is the area assigned to the point. (In general, the outward normal deriva- 
tive 94 >/ 9 n is enforced at a boundary point by applying a "force” to the unconstrained 
DOF at that point equal to G e A94>/9n. A positive force corresponds to a positive out- 
ward normal derivative.) In Eq. (19), the &2 term is analogous to a scalar spring of 
constant G e Aa 2 /ai connected between the point and ground. The 83 term is analogous 
to a scalar dashpot of constant G e Aa 3 /ai connected between the point and ground. The 
84 term is analogous to an added mass of value GgAa^ai attached to the point. 

(Here, one should probably use a consistent, rather than lumped, formulation since 
Zarda and Marcus (Ref. 30) showed that the differences between the two are not 
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insignificant for free surface flow problems.) The a$ term is a time-independent 
force given by — GgAa 5 /aj. As expected, the special case of the Neumann boundary 
condition (<|> >n = 0) corresponds to the traction-free boundary in elasticity and 
hence is a natural boundary condition. The Dirichlet condition (<j>=<() 0 ) is implement- 
ed merely by enforcing the desired value as a "displacement" boundary condition. 

SCATTERING FROM ELASTIC BODIES 

In the scattering problem (as shown in Fig. 1), a submerged elastic body is 
subjected to a plane wave time-harmonic acoustic incident loading of circular fre- 
quency a). Without loss in generality, we can assume that the waves propagate in 
the -x direction. The speed of such propagation is c, the speed of sound in the 
fluid (usually water). 



Figure 1 - The Scattering Problem 


Within the fluid region, the total fluid pressure p satisfies the wave equa- 
tion, Eq. (1). Since the free— field incident pressure p-^ is known and is given by 

p t (x,t) = p Q e 1 ^* + wt ) (20) 

(where k=to/c) it is frequently convenient to decompose the total pressure p into 
the sum of incident and scattered pressures 


P " Pi + P s 


( 21 ) 


each of which satisfies the wave equation. 

The solution of this scattering problem has been approached in various ways, 
including the T— matrix method (Ref. 16), numerical approaches using approximate 
fluid loading schemes (Ref. 15), and finite element schemes (Ref. 14). The latter 
approach models with finite elements both the structure and a portion of the infi- 
nite fluid, which is terminated with a simple radiation boundary condition to ab- 
sorb the outgoing waves as much as possible. Kalinowski's finite element approach 
(Ref. 14) uses the total fluid particle displacement as the fundamental fluid 
unknown . 
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The approach presented here is most similar to Kalinowski’s approach (of those 
mentioned above), except that we use as the fundamental fluid unknown the scattered 
component of a fluid velocity potential rather than the total fluid displacement. 
Both approaches yield symmetric matrix equations, but the use of velocity potential 
rather than displacement as fluid unknown results, in 3-D, in fluid matrices of 
one-third the order and one-third the matrix bandwidth. 

We now formulate the problem for finite element solution. The finite element 
modeling of the elastic structure results in the matrix equation 

Mu+Bu+Ku=-Ap=-A (p^ + p g ) (22) 


where u is the vector of displacement components in the structure; M, K, and B are 
the structural mass, stiffness, and damping matrices, respectively; p is the vector 
of fluid pressures at the nodes of the fluid region; and A is the area matrix which 
converts fluid pressures at interface points to structural loads. 

A finite element model of the fluid region results in a matrix equation of the 

form 


Q p s + H P s = 0 


(23) 


where p s is the vector of scattered fluid pressures at the nodes of the fluid re- 
gion, and Q and H are the fluid "inertia" and "stiffness" matrices, respectively. 
According to the analogies described in the preceding section, the same finite ele- 
ment code may be used to model both the structural and fluid regions. From Eq. 
(15), material constants assigned to the elastic elements used to model the fluid 
satisfy 


E e = 


a G„ 


= G e /c' 


(24) 


where a is given by Eq. (16). 

As is, Eq. (23) does not account for either the fluid-structure interface con- 
dition (2) or a wave— absorbing boundary condition. As in Kalinowski's work (Ref. 
14), we will use the simple plane wave absorbing condition (Ref. 31) 


3p/8n = -p/c 


(25) 


Other possibilities are discussed by, for example, Engquist and Majda (Ref. 32), 
Bayliss, Gunzburger, and Turkel (Ref. 33), and Israeli and Orszag (Ref. 34). How- 
ever, these are not of the general form of Eq. (10). Kalinowski showed (Ref. 14, 
26,35,36), and our Example 1 verifies, that the plane wave absorbing condition is 
satisfactory if the outer boundary is far enough away from the structure. 
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Both boundary conditions (2) and (25) can be handled using Eq. (19). At the 
fluid-structure interface, from Eqs. (2) and (21), 


•• •• •• 

9p s /3n = -Sp^an - pu n = p (u ni - u fl ) 


(26) 


where n is the outward unit normal from the structure (into the fluid), and (i n and 
u n i are, respectively, the total and incident outward components of fluid particle 
acceleration at the interface. Thus, using Eq. (19), we impose the condition (26) 
by applying a load to each interface fluid point given by 


F (P) - -G 0 Ap (u nl - u n ) 


(27) 


Similarly, the radiation condition (25) is enforced by applying a load to each 
fluid point on the outer boundary given by 


F (p ) = -(G A/c) p c 


(28) 


That is, a dashpot of constant G e A/c is connected between each boundary point and 
ground. 

The overall matrix system describing the coupled problem is obtained by combin- 
ing Eqs. (22), (23), (27), and (28): 


M 0 


"GgPA 0 


u 


B 0 

0 C 


u 


+ 


K A 

0 H 


u 


-Ap ± 

- G eP A %i 


(29) 


This system, which is nonsymmetric , can be symmetrized (Ref. 28) by reformulating 
the problem to use a new fluid unknown q such that 


Ps 


= q 


(30) 


If the second partition of Eq. (29) is integrated in time, if p s is replaced by 
q, and if the fluid element "shear modulus" G & is chosen as 


G e = -1/p 

the system (29) is transformed into 


(31) 


M 

0 

I 

u 


B 

A 

| 

f • 

U 



< 


| + 



i 


0 

Q 

1 

q 


a t 

C 

1 

| • 

1 <1 


K 

0 


u 


-Ap ± 

0 

H 


q ] 

= < 

Av ni 


(32) 
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where v n ^ = u n ^. This is the form of the equations which we will use to solve the 
scattering problem. 


The new variable q is, except for a multiplicative constant, the velocity po- 
tential <(> long used by fluid dynamicists, since 


P = "P<(> (33) 

Eq. (32) could also be recast in terms of <j> rather than q as the fundamental fluid 
unknown, but no particular advantage would result. In fact, the use of q rather 
than <)> has the slight practical advantage that the fluid pressure can be recovered 
directly from the finite element program as the time derivative (velocity) of the 
unknown q . 

To summarize, both structural and fluid regions are modeled with finite ele- 
ments. For the fluid region, the material constants assigned to the finite elements 
are 


G e 1/p, E e - a G e , P e = -1/pc 2 (34) 

where a is given by Eq. (16). The dashpots making up matrix C in Eq. (32) are ap- 
plied at the outer fluid boundary with damping constant -A/ pc at each grid point to 
which the area A has been assigned. At the fluid-structure interface, matrix A is 
entered using the areas (or areal direction cosines) assigned to each wet degree of 
freedom. 

The right-hand side of Eq. (32) can be simplified further since, for plane 
waves propagating in the negative x-direction at speed c, the free-field incident 
pressure and incident fluid particle velocity in the x-direction are related by 
(Ref. 37) 


Pi = -pc v xi 


(35) 


Then, as in Fig. 1, if we define 9 as the angle between the normal n and the posi- 
tive x-axis, 


v ni = v xi cos 9 (36) 


The x-component of the free-field fluid particle velocity v xi is the same at all 
points in space except for a phase angle, which may be introduced into the analysis 
by means of the time delay between two points having different x-coordinates . 

EXAMPLE 1: SCATTERING FROM INFINITE CYLINDRICAL SHELL 

The formulation derived in the preceding section will be illustrated first on 
the two-dimensional (plane strain) problem of scattering from an infinitely long 
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cylindrical steel shell (Fig. 2) of radius 
properties of steel are E = 19.5xl0 10 N/m* 
surrounding medium is seawater, for which c 


0.5 m and thickness 0.01 m. The 
v =0.28, and p g = 7700 kg/m 3 . 

- 1500 m/sec and p = 1026 kg/m 3 . 


material 

The 



PLANE 

WAVE 


Figure 2 - Scattering from Infinite Cylindrical Shell 


latera^four-node “t^s me" 1 !! ^TT th * Wlth 2 ' D ^adri- 

ing to Eq. (34). The thickness in the z-direction 3 ? or^lTi asslgned accord- 

chosen as 1 m. By symmetry only half nf £® ctio " for aj 1 elements is arbitrarily 
to be modeled. 7 ’ 7 ° f the problem (the half-plane y > 0) needs 


For a frequency of excitation of 4100 Hz fka = » m „„ , , , , 

with 128 elements spanning the 180 degrees of Mr. • ) » we modeled the region 
ary (a circle of radius 1.026 m) was 1 44 am f cumfer ® nce * The outer fluid bound- 
According to Kalinowski (Ref 141 Mol a acoast * c wavelengths away from the shell, 
absorption of outgoing ^waves! ' ' ° f fWd Should excellent 

of the surface sc^t'ere'rpJessu^s” IVi’ ” hi ' h £h °'' 8 a prison 

an analytic expression presented by lunger and Felt^r 38) el< Th" tS “ d by USi " g 
clearly excellent. Although the T,.no- 0 r-Ii?ex-t - ^ t ’ 38 '* The agreement is 

series solution (and hence another numerical apjroach^it C ° mparlson is a 
concerning the radiation boundary condition PP oach) ’ ic makea no approximation 

S ho™^“r^ b e re b«r;re::„^rbrHe b XL“ snsr both above and ^ 

SCATTERING FROM FREE— FLOODED ELASTIC SHELL WITH RIGID INNER CORE 

elast^r^ SlTj^h^ free a submerged 

the shell and also contains a rigid inner core (£“%)? “ S the OUCslde > <* 

vious^problem“except tLj h Jhe P ^ b daLntarn a id l0n v SlmlUr tC that of cha P™- 

the total velocity potential rather than the < 5 ^ 3 ^ ln the inner fluld region is 

tial as in the outer f!uid a region! r than SCattered com P onent of velocity poten- 
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90 deg 



Figure 3 - Scattered Pressures on the Surface of an Infinite 

Cylinder (ka = 8,6) 



Figure 4 - Scattering from Free-Flooded Elastic Thin Shell 

with Rigid Inner Core 


Therefore, let p ls denote the scattered pressure in the outer fluid region a 
p 9 denote the total pressure in the inner fluid region. Since the elastic shell is 
thin the same shell variables interact with both inner and outer fluid regions. 

The main difference between the treatment of the two interfaces is that the normal n 
in Fig. 4 is an inward normal for the outer fluid and an outward normal for the inner 

fluid. 
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The generalization of Eq. (29) can then be written immediately as 




' h 2 

0 

0 " 


' P2 \ 

r 0 ' 

+ 

-A 

K 

A 

4 

' ' * 

-A Pi 


0 

0 

H 1 - 


^ Pis' 

Au nl , 


(37) 


where G e has been specified as in Eq. (31), and we have assumed that the same area 
matrix A is used on both interfaces to convert pressures to forces. 


This nonsymmetric system can also be symmetrized by reformulating the equations 
in terms of two new fluid variables qj and q 2 such that 


$1 = Pis’ 9 2 = P2 (38) 


As before, the nonsymmetric terms can be moved to the damping matrix by integrating 
the first and third partitions of Eq. (37) in time: 


0 2 0 0 

0 M 0 

- 0 0 Qj 


92 ] 

0 

-a t 

0 


' 9 2 ’ 

u > + 

-A 

B 

A 

< 

A 

*1 ' 

. 0 

a t 

C _ 


< j 



' h 2 

0 

0 ’ 

( 92 1 

( 0 ■ 

+ 

0 

K 

0 

1 “ 1 

" ■ “API 


. 0 

0 

Hi - 

l 91 \ 

l Av ni , 


(39) 


where, for both fluid regions, the material constants assigned to the finite ele- 
ments are given by Eq. (34). 
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EXAMPLE 2: SCATTERING FROM INFINITE CYLINDRICAL SHELL ENCLOSING A FLUID 

BOUNDED BY A CONCENTRIC RIGID CYLINDER 

The formulation derived in the preceding section will be illustrated on the 
two-dimensional (plane strain) problem of scattering from an infinitely long cylin- 
drical steel shell which contains on the inside both fluid and a concentric rigid 
cylinder (Fig. 5). The rigid inner core has a radius of 0.254 m. The shell is 
identical to that used in Example 1. Seawater floods the region between the shell 
and the rigid inner core. 


STEEL 



Figure 5 - Scattering from Infinite Cylindrical Shell Enclosing 

a Fluid Bounded by a Concentric Rigid Cylinder 


A finite element model was prepared for this problem for excitation at 2100 Hz. 
(ka = 4.4). Since this frequency is lower than that used in Example 1, a coarser 
mesh (having 96 elements spanning the 180 degrees of circumference) was used. Here, 
the outer fluid boundary was a circle located 1.32 m away from the shell, a distance 
equal to 1.85 wavelengths of the incident free-field acoustic wave. 

For this problem, no analytic results were readily available for comparison, so 
instead a comparison was made to the corresponding problem with the inside of the 
shell evacuated (as in Example 1, but at a different frequency). See Fig. 6. 
Clearly, the presence of the contained fluid and the inner core affect the solution 
significantly. 

Although the two cylinders in this problem were chosen to be concentric, the 
formulation derived in the preceding section is clearly general enough to handle 
arbitrary geometry. In addition, the formulation applies to both 2— D and 3— D prob 
lems. 


RADIATED PRESSURES 

A finite element model of an exterior fluid-structure interaction problem is 
capable of predicting only near-field acoustic pressures because of the approximate 
nature of the radiation boundary condition on the fluid. However, given fluid pres- 
sures and normal velocities on the fluid-structure interface, the fluid pressure at 
any point in the exterior field can be calculated by a numerical quadrature. 
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— elastic cylinder empty 

— ELASTIC CYLINDER ENCLOSING 
CONCENTRIC RIGID CYLINDER 
WITH FLUID BETWEEN 


90 deg 



Figure 6 - Comparison of Surface Pressures Scattered by an Evacuated 

Cylinder with Pressures Scattered by a Cylinder Containing a 
Fluid and a Concentric Rigid Cylinder (ka = 4.4) 


In Fig. 7, let z_ be the position vector to an exterior field point P, and z = 
|zj . Let x be the position vector to a point on the fluid-structure interface 
(with x = Tx|), let _r = z_ - x (with r = |rj), and let £ be the unit outward 
normal at the location x. The pressure at z is (Ref. 39) 


p(z) = - / g q(x)(e -ikr 74irr )dS + J g p(x) 3/9n(e ikr /4irr)dS 


(40) 


where 


q = 3p/3n - -ioopv n 

(41) 

a 

3 

n 

(42) 


and it is assumed that the harmonic time-dependence of the variables is exp(ioot) 
rather than exp(-ia)t). 

Since 


3/3n(e ^ kr /4irr) = (e ^ kr /4irr)(ik + r 1) cos 3 


(43) 



p 



Figure 7 - Notation for Far-Field Radiated Pressure Calculations 


it follows that 


p(z^ = / s [io>pv n (x) + (ik + r *) p(x) cos g] (e _i ^ r /4nr)dS 


(44) 


This expression is valid for any point in the exterior field. 

Eq. (44) can be simplified if only far-field locations are of interest. As 
|z| -*■ 00 , ik + r“l > ik, and, from the law of cosines, r z - x cos a. Thus, at 
far-field locations. 


p(z) = (ik e~^^ z /4irz) |g [pcv n (3c) + p(x) cos g]e^ x cos 01 dS (45) 


where 


cos g = ( zj |£| ) • n 


(46) 


The numerical integrations in Eqs. (44) and (45) require, for each wet struc- 
tural point, the knowledge of the location (coordinates), normal, area, pressure, 
and outward normal velocity. All these quantities can be obtained directly from 
NASTRAN using the 0UTPUT2 utility module. The grid point coordinates and areas 
(or areal direction cosines) can be obtained from an abbreviated static analysis 
in which a unit outward pressure load is applied to the structure. The unit 
normal at a point is then the unit vector parallel to the area vector. The 
pressures and velocities are obtained from the frequency response analysis. 
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